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Hi Absorption Toward the Nucleus of the Radio Galaxy 
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ABSTRACT 

We present sensitive, high resolution Very Long Baseline Array (VLBA) 
observations of the central 0.3" of PKS 2322—123 at 1.3 and 5 GHz. These observations 
reveal straight and symmetric jets emerging from both sides of an inverted spectrum 
core. The 21 cm line of atomic hydrogen is detected in absorption against the core 
and eastern jet with substantial opacities, but is not seen towards the equally strong 
western jet. Both a narrow (110 km s" 1 FWHM) and a very broad (735 km s 
FWHM) line are seen, although the very broad line is seen only against the core. Both 
lines are redshifted (~220 ± 100 km s _1 ) with respect to the systemic velocity. The 
most likely explanation for the observed Hi kinematics are an atomic torus centered 
on the nucleus with considerable turbulence and inward streaming motions. The scale 
height of this torus is less than 20 pc. Though rare in flux limited samples of compact 
radio sources, symmetric parsec-scale structure appears nearly ubiquitous among radio 
galaxies with Hi absorption, probably because they are viewed more nearly edge-on 
through the torus. 

Subject headings: cooling flows — galaxies: clusters: individual (A 2597) — galaxies: 
jets — radio continuum:galaxies — galaxies: elliptical and lenticular, cD: — radio lines: 
galaxies — galaxiesdndividual (PKS 2322— 123) 
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1. Introduction 

Hi gas has now been detected towards several radio galaxies embedded in X-ray cooling 
flow clusters: 3C84 in A 426 (Crane et al. 1982; Jaffe 1990); Hydra A in A 780 (Taylor 1996; 
Dwarakanath, Owen, & van Gorkom 1995), and recently PKS 2322-123 in A 2597 (O'Dea, Baum 
& Gallimore 1994). In an even larger number of such radio galaxies Hi has not been detected in 
spite of fairly intensive searches (e.g. Virgo A, and 3C 338 in A 2199 - Dwarkanath, van Gorkom, 
& Owen 1994; Jaffe 1991; O'Dea, Gallimore, & Baum 1995). Besides a natural curiosity about 
these unusual systems, additional motivation to search for Hi has come from the discovery of 
excess low-energy X-ray absorption (White et al. 1991) which has been suggested to be due to a 
population of cold Hi clouds in the cluster with a column density of Nh ~ 10 21 cm -2 and covering 
factor near unity. To date this population has not been seen and the Hi gas that is detected seems 
to be closely associated with the active nucleus. In Hydra A the Hi gas appears to be distributed 
in a disk having a rotation axis closely aligned with the jet axis and with a height of just 30 pc. 

O'Dea, Baum & Gallimore (1994) detected two Hi components in PKS 2322—123 using the 
VLA in its A configuration. One component is a narrower (220 km s _1 wide) feature seen only 
towards the center of the source (the location of the nucleus). The other component is broader 
(410 km s _1 wide) and spatially extended against the entire radio source. O'Dea et al. suggest 
that the "broad" component is associated with a bright Ha nebula. 

A detailed study of the optical spectrum from the central ~kiloparsec of A 2597 by Voit & 
Donahue (1997) allowed them to measure the density, temperature and metal abundances of the 
line-emitting gas and conclude that the most likely explanation for the brilliant ionizing nebula in 
A 2597 is a population of hot stars. As a byproduct of their study they determine the redshift 
of PKS 2322-123 to be 0.0821 ± 0.0001. HST observations by Koekemoer et al. (1998) reveal a 
series of blue knots and clumps that are most likely regions of recent star formation. They also 
find emission line filaments at the edges of the radio lobes, which argues for a direct interaction 
between the radio source and its environment. Falcke et al. (1998) recently detected molecular 
hydrogen in the central kiloparsecs of PKS 2322—123. They suggest that the luminous H2 emission 
correlates with the presence of a strong radio source, luminous Ha emission and a cooling flow of 
>100 M y- 1 . 

We assume Hq = 75 km s^ 1 Mpc" 1 and qo = throughout. 
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2. VLBA Observations 

The VLBA observations were carried out at 1.312 and 4.991 GHz with the 10-element VLBA 
and 27-element VLA of the NRAOQ in a single 10.5-hour observing session on 1996 December 
7. The VLA was operated in "phased-array" mode to synthesize a single aperture with a large 
collecting area. The calibrator J2246-1206 was used to phase up the VLA every 25 minutes. The 
net integration time on PKS 2322— 123 was 323 minutes at 1.3 GHz and 95 minutes at 5 GHz 
spread over multiple 19 minute snapshots to improve u,v coverage. Both right- and left-circular 
polarizations were observed with 2-bit sampling. The 1.3 and 5 GHz data were correlated in 
separate passes so as to provide high spectral resolution (1024 channels per 8 MHz wide IF band) 
without the RL and LR cross-correlation products at 1.3 GHz, and low spectral resolution (16 
channels per IF band) with full cross-correlation products at 5 GHz. 

Standard a priori flagging, amplitude calibration, fringe fitting, bandpass calibration (using 
3C84), and frequency averaging procedures were followed in AIPS. Global fringe fitting was 
performed on J2246-1206 and the resulting delays and rates were transfered to PKS 2322— 123. 
Furthermore, after imaging, J2246-1206 was used in conjunction with CALIB to obtain time- 
dependent amplitude corrections which were also applied to A 2597. Subsequently, A 2597 was 
phase-only self-calibrated with a 1 minute solution interval. No amplitude self-calibration was 
performed on A 2597 itself. All manual editing, imaging, deconvolution, and self-calibration were 
performed using DlFMAP (Shepherd, Pearson, & Taylor 1994, 1995). Polarization calibration was 
performed on the 5 GHz data using J2246-1206 to determine the instrumental leakage terms, and 
the absolute polarization angle calibration was determined from a single 12 minute observation of 
3C286 which has known polarized structure (Cotton et al. 1997). 

3. Continuum Results 

In Fig. 1 we show the naturally weighted 1.3 and 5 GHz VLBA images. Both images display 
a clearly symmetric structure along a position angle of ~70°. On kiloparsec scales a bright hot 
spot 0.5" south-west of the nucleus (Sarazin et al. 1995) suggests some bending of the jet to a 
position angle of ~50°. At the hot spot the south-west radio jet appears to be sharply deflected 
to the south, giving the large scale structure an overall "C" shape. 

In the parsec scale jet at 1.3 GHz there are 5 components extended nearly in a line over 63 
mas (89 pc) with the strongest component at the center. These components are labeled in Fig. 1, 
and parameterized by Gaussian modelfitting in Table 1. The central component is considerably 
stronger at 5 GHz and has a spectral index, a, of 0.6 ±0.1 where S u cx v a . We identify this 
component as the core based on its inverted spectrum, compactness, strength, and the symmetric 
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morphology of the source. No significant polarized flux density was measured at 5 GHz above the 
noise floor of 40 /uJy/beam. This indicates a fractional polarization of <0.1% for the nucleus, and 
<10% for components E2 and W2. 



TABLE 1 

Component Flux Densities and Spectra 



Component 


r 


9 


Sl.3GHz 




#5GHz 


oti. 


3-5 




(mas) 


(°) 


(mJy) 




(mJy) 






El 


33 


64 


2.8 ± 0.5 




< 0.2 


< - 


1.9 


E2 


16 


67 


6.7 ± 0.6 


4.8 


± 0.3 


-0.25 ± 


0.1 


Core 








16.6 ± 0.9 


36.9 


± 1.9 


+0.60 ± 


0.1 


W2 


16 


-112 


7.6 ± 0.6 


4.1 


± 0.3 


-0.47 ± 


0.1 


Wl 


30 


-115 


2.4 ± 0.5 




< 0.2 


< - 


■1.9 



A spectral index map of the nucleus of PKS 2322—123, made by tapering the 5 GHz image to 
the 1.3 GHz resolution, is shown in Fig. 2. The two images were aligned using the relatively steep 
spectrum component W2, rather than the core, because if the core component is the location in 
the jet where the synchrotron self absorption optical depth is unity, then its position will vary 
with frequency (Blandford & Konigl 1979), such that higher frequencies look "deeper" down the 
jet towards the true center of activity. Also, a gradient in the free-free opacity towards the core 
may result in a shift in the apparent core location as a function of frequency. Given the low 5 
GHz brightness temperatures for E2, C and W2 of 1.4 xlO 4 , ~ 10 7 , and 1.2 xlO 4 K respectively, 
synchrotron self absorption is unlikely to be important. Higher frequency VLBI observations are 
needed to determine the intrinsic flux density and spectral index before the free-free optical depth, 
Tf f, can be measured, but a rough estimate on Tff can be obtained by extrapolating the 5 GHz 
flux density to 1.3 GHz assuming a spectral index of —1. This gives Tff of 1, 2.1 and 0.7 for E2, 
C and W2 respectively. These optical depths require emission measures of 7 xlO 7 , 14 xlO 7 , and 
5 xlO 7 cm~ 3 pc for E2, C and W2 assuming a gas temperature of 10 4 K (see Levinson, Laor, & 
Vermeulen 1995). 



4. Hi Absorption Results 

In Fig. 3 we show the 1.3 GHz continuum emission and the spectra of each of the 5 
components. Both a very broad (735 km s" 1 ) and a narrow (106 km s _1 ) line are seen against the 
core component. Only the narrow line is seen against component E2, and neither line is detected 
towards W2. Measurements of the lines are given in Table 2 and shown graphically in Fig. 4. 
Both of these lines are significantly redshifted with respect to the velocity centroid of the optical 
emission lines from the nucleus. Voit Sz Donahue (1997) have obtained a high signal-to- noise 
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TABLE 2 
HI Column Densities 



Component 


5l.3GHz 


S'line 


V 


AV 


T 






N H 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 






(7) 


El - narrow 


2.8 ± 0.5 


< 1.5 mJy 


- 


- 


< 0.8 






- 


E2 - narrow 


6.7 ± 0.6 


3.1 ± 0.3 


24882 ± 6 


133 ± 14 


0.63 


12.2 


X 


10 23 


Core - narrow 


16.6 ± 0.9 


2.3 ± 0.6 


24880 ± 12 


106 ± 33 


0.15 


2.3 


X 


1Q 23 


W2 - narrow 


7.6 ± 0.6 


< 2.0 mJy 






< 0.3 








Wl - narrow 


2.4 ± 0.5 


< 1.5 mJy 






< 1.0 








El - broad 


2.8 ± 0.5 


< 1.0 mJy 






< 0.5 








E2 - broad 


6.7 ± 0.6 


< 1.0 mJy 






< 0.2 








Core - broad 


16.6 ± 0.9 


3.8 ± 0.3 


24905 ± 19 


735 ± 53 


0.26 


27.8 


X 


10 23 


W2 - broad 


7.6 ± 0.6 


< 1.0 mJy 






< 0.2 








Wl - broad 


2.4 ± 0.5 


< 1.0 mJy 






< 0.5 









Notes to Table 2 

Col.(l). — Component name and type of Hi line detected (or upper limit given for). Col. (2). — 
Continuum flux density at 1.3 GHz in mJy. Col. (3). — Depth of the line in mJy or 3a upper limit. 
Col. (4). — Central velocity in km s -1 . Col. (5). — FWHM in velocity in km s _1 . Col. (6). — Optical 
depth or 3a upper limit. Col. (7). — The column density in units of cm~ 2 calculated assuming a 
covering factor of 1 and a spin temperature of 8000 K (Conway & Blanco 1995). 

spectrum from which they derive an emission line redshift of 24613 ± 60 km s _1 . The presence 
of large scale asymmetries in the emission line gas kinematics, however, may indicate that the gas 
is substantially offset from the systemic velocity. Stellar absorption features may better trace a 
more relaxed inner potential. Taking advantage of prominent Ca II H, K, and G band absorption 
features we performed a cross correlation analysis of the Voit & Donahue spectrum with a range of 
galaxy template spectra. Correlations with different spectra gave consistent results with a mean 
systemic velocity of 24673 ± 20 km s _1 . Given the peculiar nature of PKS 2322— 123 (elliptical 
with strong emission lines), which is not well matched to any of the templates, and the mismatch 
in spectral resolution, a more realistic error is ±100 km s _1 . The mean offset of the Hi absorption 
features from the systemic velocity is thus 220 ± 100 km s -1 . 

The broad (412 km s" 1 FWHM) Hi line centered at 24604 ± 17 km s" 1 discovered by O'Dea 
et al. (1994) is not seen towards any component on the parsec scale. This is not surprising since 
it was spatially resolved, and appeared strongest towards the NE lobe of PKS 2322— 123. In all 
likelihood the continuum and Hi emission from the NE lobe has been resolved out. The narrow 
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(221 km s 1 FWHM) spatially unresolved line seen by O'Dea et al. at 24886 ± 5 km s 1 towards 
the nucleus is detected and spatially resolved by our observations. 

The very broad (735 km s _1 FWHM) line was not detected by O'Dea et al.. Its depth of 3.8 
mJy should have been detectable by them even when diluted by unabsorbed emission within the 
beam, but given their poorer velocity coverage (1304 km s _1 after editing) it is possible that this 
very broad line was lost in the continuum subtraction. 

5. Discussion 

5.1. Relation of Symmetric Structure to Hi Absorption 

The parsec scale structure of PKS 2322— 123 appears extremely symmetric at both 1.3 and 5 
GHz. If the jet components are intrinsically similar and were ejected from the core simultaneously, 
then cither (1) the radio source must lie within a few degrees of the plane of the sky or (2) the bulk 
motion of the jets must be nonrelativistic. VLBI observations of two other FR-I radio galaxies 
embedded in X-ray cooling flow clusters (Hydra A - Taylor 1996; and 3C 338 - Feretti et al. 
1993) have revealed remarkably symmetric structures on the parsec scale whereas over 95% of all 
AGN in complete flux limited samples have one-sided core+jet morphologies (Taylor et al. 1996). 
If even low power radio galaxies have jets that start out relativistic on parsec scales (Giovannini 
et al. 1998), then the jets of all 3 sources must be oriented with a few degrees of the plane of the 
sky. The chance of 3 out of 13 radio galaxies embedded in cooling flow clusters (see Taylor, Barton 
& Ge 1994) being within 5 degrees of the plane of the sky is ~3%. Alternatively, it has been 
suggested (e.g., Soker & Sarazin 1988) that cooling flows may decelerate the radio jets, although 
such disruption has been considered to occur at the sonic radius (where the inflow velocity exceeds 
the sound speed) somewhere between 0.1-10 kpc. Unless the cooling flows have an affect on the 
jet closer to the nucleus than previously thought it is difficult to understand why the parsec scale 
jets in Hydra A, 3C 338, and PKS 2322— 123, would have such low bulk velocities. 

Examination of the parsec scale structure of sources with Hi detected in absorption (Hydra A 
- Taylor 1996; 1413+135 - Perlman et al. 1996; N3894 - Peck & Taylor 1998; 1946+708 - Peck, 
Taylor, h Conway 1999) indicates that a symmetric parsec-scale structure correlates strongly with 
the presence of Hi, although there are counter-examples (e.g. Mrk231 - Carilli, Wrobel & Ulvestad 
1998). 3C84 could be an intermediate case - observations by Vermeulen et al. (1994) suggest an 
inclination of ~45°. The most likely explanation for this correlation is that symmetric sources are 
oriented close to the plane of the sky such that their appearance is not made one-sided by Doppler 
boosting effects. This then provides a favorable viewing angle of the radio source through a disk 
or torus of atomic gas that is perpendicular to the jet axis. This picture is consistent with unified 
schemes, and our limited understanding of how the central engine in AGN might be fed. 

In the inner regions of the disk a large fraction of the gas must be ionized by the central 
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engine. This will result in free-free absorption of the radio continuum at frequencies below ~5 
GHz as seen in 3C 84 (Walker et al. 1998), Hydra A (Taylor 1996), 1946+708 (Peck et al. 1999), 
and PKS 2322— 123. Another result of the dense ionized gas, if magnetized, could be extremely 
high Faraday rotation measures (RMs). Owing to the lack of any polarized flux in these systems, 
it has not been possible to directly measure the RMs in any source with Hi absorption, but the 
presence of very high RMs could explain why no polarized flux is detected from these sources. 

5.2. Location of the Atomic Gas 

The sharp spatial variations in optical depth on scales of order 10 pc (Fig. 4), and the 
detection of the very broad line only against the core, argue that the atomic gas we are seeing is 
located in the central 10s of parsecs of PKS 2322— 123, and is not just a chance encounter with a 
cloud in the galaxy. If the gas is distributed in a disk centered on the core then it must be fairly 
thin (<20 pc) so as to cover the core and E2, but not W2. This situation is somewhat different 
from the radio galaxy 1946+708 (Peck et al. 1999) where Hi gas is seen across the entire 60 
parsec extent of the radio source. However, in 1946+708 as in PKS 2322— 123, a broad component 
(~300 km s -1 wide) is seen only within ~10 pc of the core. The 735 km s _1 FWHM line in 
PKS 2322— 123 is, to our knowledge, the broadest Hi absorption line, and is considerably wider 
than the optical emission lines from the nucleus of PKS 2322—123 which have a velocity dispersion 
of 270 km s _1 (Voit & Donahue 1997) albeit from much lower spatial resolution observations. 

The fact that both lines in PKS 2322-123 are redshifted by ~220 ± 100 km s^ 1 may imply 
that the gas is in-falling. This is similar to the trend found by van Gorkom et al. (1989) for 6 of 
8 Hi absorption systems in a sample of nearby ellipticals. If this gas is in a disk, then the disk 
is not in an orderly Keplerian rotation about the nucleus, but has significant inward streaming 
motions. From the very broad linewidths, much broader than the thermal linewidths for any 
plausible conditions, we can also infer considerable velocity structure in the disk, probably due to 
turbulence. Alternatively, given the velocity offset of the very broad line from the absorption line 
redshift of 232 ± 100 km/s, it is possible that the very broad line is at the systemic velocity. In 
this case it would be natural to explain the linewidth by rotation. Unfortunately, we don't know 
the radial extent of the disk since that is perpendicular to the radio axis. Assuming a radius of 10 
pc, the mass required to produce the observed linewidth is 2 xlO 8 MQ. 

GBT thanks STScI for hospitality during a visit in which this paper was written. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by 
the Jet Propulsion Laboratory, Caltech, under contract with NASA. 
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Fig. 1. — The naturally weighted images of PKS 2322—123 at 1.3 and 5 GHz. The restoring beam 
has dimensions 16.3 x 6.7 mas in position angle —1.3° at 1.3 GHz and 4.63 x 1.86 mas in position 
angle —5.1° at 5 GHz. Contours are at factor 2 intervals and start at 0.5 and 0.2 mJy/beam for 
the 1.3 and 5 GHz images respectively. 
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Fig. 2. — A spectral index map of PKS 2322-123 between 1.3 and 5 GHz with 1.3 GHz contours 
overlaid. The 5 GHz image has been heavily tapered and restored to the same resolution as the 1.3 
GHz image. The images have been registered using component W2. 




100 -100 

Relative R.A. (mas) 

Fig. 3. — The 1.3 GHz continuum image from Fig. 1 with the spectra for each of the 5 components 
overlaid. The velocity resolution is 7.1 km s _1 . The arrow at 24673 km s _1 indicates the systemic 
velocity derived from the optical redshift as discussed in the text and the width of the arrowhead 
represents the 1 a uncertainty of 100 km s _1 . 



